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Tj- ' Abstract 

O ' 

We study the heavy charge potential in the Coulomb phase of pure gauge compact U(l) 
theory on the lattice. We calculate the static potential Vw(T,R) from Wilson loops on 
a 16 3 x 32 lattice and compare with the predictions of lattice perturbation theory We 
investigate finite size effects and, in particular, the importance of non-Coulomb contributions 
to the potential. We also comment on the existence of a maximal coupling in the Coulomb 
phase of pure gauge U(l) theory. 

^ ■ 

H : 1 Introduction 



In 1974 Wilson proposed the compact lattice formulation of pure gauge U(l) theory The 
action of this model is 

M^) = /?EE(i-cos^), P = \, (i) 



X fl>U 



where g 2 is the bare coupling constant, and the link variables are U Xfl = exp(i9 x ^), 
Oxfi e (-7r,7r]. The plaquette angles are given by 6 XtfMU = 9 X)IX + x +p,,v - x +o^ ~ Q x ,v 
This action makes up the pure gauge part of the full QED action Sqed, which is supposed 
to be compact if we consider QED as arising from a subgroup of a non-abelian (e.g., grand 
unified) gauge theory 0. 

* Work supported by the Deutsche Forschungsgemeinschaft under research grant Mu 932/1-4 and EEC- 
contract CHRX-CT92-0051 and by the US DOE under grants # DE-FG05-85ER250000 and # DE-FG05- 
96ER40979. 

^present address: BLTPh, JINR, Dubna, Russia 
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The average of any gauge-invariant functional 0(U) is 

(0(U)) = Z- 1 [HdU x » 0(U) ■ e~ s ^ , dU Xfl = ^ , (2) 

where Z is the partition function, defined formally by (1) = 1. At small enough /3's the 
strong coupling expansion shows an area-law behaviour of the Wilson loops, while at large 
/3's a deconfined phase exists. The two phases are separated by a phase transition at some 
'critical' value j3* whose existence was assumed in The weak coupling - deconfined - 
phase is expected to be a Coulomb phase, i.e., a phase with massless noninteracting vector 
bosons (photons) and Coulomb-like interactions between static charges. 

There are scarcely any doubts about the existence of the Coulomb phase in the weak 
coupling region for Wilson's QED, though - to our knowledge - there is no rigorous proof. 
It is worthwhile to mention that for the Villain approximation 0] such a proof is available 
[f|. [H . However, a detailed study has shown that the Villain action is quantitatively a rather 
bad approximation to the Wilson action in the weak coupling region 0. 

A perturbative expansion suggests for the lattice potential V latt (R) the expression 

oo 

V lat \R) = Y, 9 2n ■ V n (R) , (3) 

71=1 

where up to an additive constant V\ is the lattice analog of the 

continuum Coulomb potential ~ 1/AirR (for an explicit expression see fl7|) below). One- 
loop corrections (~ g 4 ) do not change the functional dependence of the potential but at the 
two-loop level (~ g 6 ) non-Coulomb-like contributions appear (see below). The analytical 
and numerical study of these contributions, as well as of the finite volume behavior of the 
potential, constitutes the aim of the present work. 

The behavior of the heavy charge potential in the weak coupling phase in pure U(l) 
gauge theory was the subject of several numerical studies (see, e.g. ff\, |8|, 3, |I0[|). In all 
cases consistency with a Coulomb-like behavior was reported. However, we feel that a more 
elaborate and systematic study is necessary. In particular, Monte Carlo simulations give 
precise enough measurements of the potential, that finite size effects and finite lattice spacing 
effects, i.e., deviations of the lattice Coulomb potential from the continuum 1/R behavior 
should not be neglected. In addition, finite T effects in the extraction of the potential from 
Wilson loops have to be taken into account. We shall find, on the other hand, that the 
non-Coulomb-like contributions to the potential are negligible. 

We consider this systematic study of the Coulomb phase in pure gauge U(l) theory a 
necessary step before attempting an investigation of lattice QED with fermions with the aim 
of obtaining a continuum limit that reproduces weak coupling continuum QED. 

The second section is devoted to the two- loop perturbative study of the potential, Vp(R), 
as extracted from Polyakov loop correlations. Special emphasis has been put on the study 
of the non-Coulomb-like contributions that appear at two-loop order, i.e., at order 0(g 6 ). 
In the third section we discuss the potential, Vw(T,R), as extracted from Wilson loops. 
We present the results of Monte Carlo simulations and their comparison with perturbation 
theory. Indeed, perturbation theory will play a major role in our extraction of a renormalized 
coupling, g\, from the numerical simulations. In section ^ we address the question of 
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Figure 1: Different diagrams contributing to the potential Vp(R). 



the existence of a maximal coupling in pure gauge U(l). The last section is reserved for 
conclusions and discussions. 



2 The potential from Polyakov loop correlations at 

0(f 



One way to define the heavy charge potential V = Vp(R) is 



V P (R) = -±--\nT P (R) , (4) 



where we consider a finite lattice of size Ng x N4 and Tp(R) is the Polyakov loop correlator 

I f e wE./^ 

Z J A 

if the currents correspond to a static heavy charge-anticharge pair 



T P (R) = ^ / e i9 ^ J *» A *» ■ e~ SwiA) , (5) 

Z J A 



rP( T \ _ j % - at < X 4 < iV 4 - 1 , fi v 

Ji {X) ~ \ otherwise ' 1 J 

and Jf (x) = 0, i = 1,2,3. 

We calculated the potential Vp(R) perturbatively up to order 0(g 6 ). Graphically the 
different contributions are shown in Figure |l|. 

In the tree approximation (Figure 0a) one obtains the lattice Coulomb potential 

9 Vcoul \R) - TPs 2^, ^ ' 
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where /Q = 2sin^-. In the infinite volume limit, N s — > oo , this potential becomes 

9 Vcoui( R ) = y_^p ^ . ^ — > 00 » (8) 

The continuum analog of this expression is 

g 2 V cont {R) = --rs + Const . (9) 

The lattice Coulomb potential defined in eq. (0) becomes close to the continuum expression 
in eq. (§) when 1 -C -R <C |iV s for sufficiently large lattice size N s . For lattice sizes that are 
typically used in numerical simulations (i.e. N s ~ 16 -j- 32) the lattice Coulomb potential in 
eq. (0) differs considerably from the continuum expression in eq. (Q). 

The 0(g A ) contribution (Figure |TJb) as well as 0(g 6 ) contributions shown in Figures [Tjc 
to [T|e result in a renormalization of the coupling without changing the form of the potential. 
However, two diagrams at order ~ 0(g 6 ) contain non-Coulomb-like contributions. One of 
them is the four-prong-spider graph (Figure |l|f ) which gives the contribution to the potential 



vU \ r ) = -^Ew^n^v^n^- ( io ) 

6I ^s {p-W} i=l 1 j=l i=l/C 2 (pW) 

Another non-Coulomb-like contribution comes from the two-loop bubble diagram shown in 
Figure |I|g 



V <*)(% = ^V Coul (R) + V^ oul (R) , (11) 

s p [1^ ) i,j=l 

where 

TiM = j^y 2 E D ij(q) D ij(q') D ij(p -q-q'), A;(<z) = «i±M^ (g) . ( i 3) 

Therefore, up to order 0(g 6 ) the static charge potential Vp(R) is 

V P (R) = g\_ loov ■ Vcoui(R) + g 6 ■ V nCo ui(R) , (14) 

with 

V nCo ui{R) = W\R) + V$ ovI (R) , (15) 

and where 

gli oop = g 2 ^ + -/ + f/) ■ 

To get an idea of the possible importance of the non-Coulomb term, we computed it 
by numerically performing the necessary lattice momentum sums on several lattices up to 
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Figure 2: Comparison of the non-Coulomb part V n couh ec L- (|T5|) , of the potential to the 
(lattice) Coulomb part Vcw, eq. (0), on a 16 4 lattice. 



size 16 4 .| As an example we compare V n c ou i with Vcoui in Figure |[ We see that V n coui is 
almost three orders of magnitude smaller than Vcoui, a difference that will even be enhanced 
when including the coupling constants for weak coupling. Furthermore V n c ou i approaches 
its asymptotic value much faster than 1/R. Since the long distance behavior is governed 
by the small momentum region, we can estimate it in the infinite volume limit from power 
counting in the small momentum region of the integrals. For the two-loop bubble diagram, 
eq. (|12|) , since the zero-momentum part corresponding to 7^(0), which contributes to Vcouh 
was split off, we expect a 1/R 3 approach to a constant. For the four-prong-spider graph, 
eq. (|10D, we expect, apart from a constant, a 1/R 5 fall-off at large distance. Log-log plots 
of the numerically computed non-Coulomb contributions on finite lattices confirm these 
expectations. 

From the large distance behavior in the infinite volume limit we can also estimate the 
finite size effects. Periodic boundary conditions can be mimicked by mirror charges, e.g. for 
an on-axis distance the infinite volume 1/R n becomes in a finite volume with box size iV s , 



1 1 

~RF + (N s - R) n 



1 1 

1 hn 

R n 



R 1 / r^ N 



(17) 



3 To perform the necessary lattice sums took a couple of weeks CPU time on an IBM RS6000 workstation 
for the largest lattice! 
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Therefore we expect the leading finite volume corrections to affect the constant at 0(1/N S ) 
for Vcoui and only at 0(l/iV|) for V n coui- The finite size effects in V n c ou i appear to be much 
smaller than those in Vcoui- Our numerical results are in agreement with this expectation. 



3 The potential and the renormalized coupling from 
Wilson loops 

The potential can also be obtained from Wilson loops: V = Vw(T, R). It is defined as follows 

V W (T;R) = In W ( T > R \ (18) 
W(T+1,R) 

— * — * 

where W(T, R) is the Wilson loop with 'time' extension T and space extension R. 

The Wilson loop can be on-axis or off-axis, and the space-like parts of the loop can 
include the contribution of many different contours. In our calculations we have chosen 
planar loops and two types of non-planar contours: the 'plane-diagonal' contour with space- 
like part in the plane (£2,^3) at fixed X\ (say, X\ = 0) connecting points x<i = X3 = and 
X2 = £3 = Ro, and the 'space-diagonal' contour in the 3d space (x\, x%, 23) connecting points 
x\ = X2 = X3 = and x\ = 22 = ^3 = -Ro- For the off-axis loops we average over all paths 
that follow the straight line between the endpoints as closely as possible. For example, for 
the loop through the origin and x 2 = x 3 = R the paths considered go through all points 
with x 2 = ^3 i n between. We average over all combinations of first taking a step in the 
2-direction followed by a step in the 3-direction and vice versa. 

At lowest order in perturbation theory the potential Vw(T; R) becomes 

V W (T;R) =w(T,R) -w(T+l,R) , (19) 

with w(T,R) the 0(g 2 ) part of the Wilson loop, which in turn, in Feynman gauge, can be 
split as 

w{T, R) = w s (T, R) + w t {T, R) , (20) 

where w t and w s are contribution of the time-time ('electric') and space-space ('magnetic') 
parts, respectively. The time-time contribution is 

^ *> - ~ A P - -Pfi) • • G M , (2!) 

where G (p) = 1//C 2 . It is easy to show that in the limit N4 — > 00 , T^oo, T <^ jV 4 the 
'time'-like ('electric') part of the potential gives the lattice Coulomb potential as obtained 
from Polyakov loop correlations and given in eq. ([?[). 

For the planar Wilson loops in the plane (23, £4) the space-space contribution w s is 

w s = -jf^ E(l - oos^ ■ l ~ C0 ^ R . Go{p) . (22) 

N * N I p^o 1 - COS P3 
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For the 'plane-diagonal' contour, in the plane (x 2 , 23) at fixed x\ connecting points 
x 2 = £3 = with x 2 = £3 = -Ro the space-space contribution is 



9 v-/o, , \ 1 - cosR (p 2 +p 3 ) . 

2N A Nf ^ 1 - cos(p 2 + Pa) 

and finally for the 'space-diagonal' contour connecting points X\ = x 2 = £3 = with 
X\ = x 2 = x 3 = Rq the space-space contribution is 



2 9 2 v- a(-\ 1 ~ cosR o(Pi +P2 +Ps) , 

2^ A{P) ■ — ; / : : r~ • (1 - cos p^T) ■ G {p) 



N i N s n^O 1 ~ COS (Pl +P2+ P3; 



A 



1 ' , P2+P3 , P2 - P3 1 2 
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2 cos h cos ; 



+ jpermutations of Pi,P2,Ps} • (24) 



The contribution from w s to Vw{T;R) vanishes in the limit T — > 00 and in this limit 
the potentials extracted from Polyakov line correlations and from Wilson loops agree. At 
finite T, however, w s gives a contribution « cR/T 2 to VV(T; i?). At a finite T the potential 
Vw(T]R) hence appears to have a small confining contribution. Vw(T;R) therefore needs 
to be either carefully extrapolated to infinite T, or the finite T effect has to be taken into 
account in the analysis, as we shall do in this paper. 

Just as was the case for the potential obtained from Polyakov loop correlations, at one- 
loop level the only effect is a renormalization of the coupling g 2 — > g 2 (l + \g 2 )- At two 
loops in the perturbative expansion we only considered, for simplicity, planar Wilson loops. 
We find the same structure as for the potential from Polyakov loops, eq.s (|14"D and (0). In 



particular, the two-loop renormalization of the coupling is identical, as given in eq. (jig). 

Our Monte Carlo data were produced on a 16 3 x 32 lattice for 11 values of (3 shown in 
Table |]. We measured planar and nonplanar Wilson loops W(T, R) with T < T max = 15 
and R < R ma x = 8v3. To decrease the statistical noise we have used the Parisi-Petronzio- 
Rapuano trick [llj for all time-like links. Because of this we do not have a valid measurement 
of the potential at distance R — 1, and this distance is therefore excluded from all further 
considerations. 

In the Coulomb phase the potential obtained from Monte Carlo simulations is expected 
to be very similar to the perturbative lattice Coulomb potential. Since the latter can have 
sizable finite size and finite T effects, as discussed earlier, we decided to use as a fit formula 
for the numerical potential Vw(T;R) the lowest order perturbative expression defined in 
eq.'s (|l9|) to fl24]) , with the only fit parameter being the renormalized coupling constant g\. 
The success of such fits will indicate that the non-Coulomb contributions are small, like in 
perturbation theory, and that non-perturbative effects just go into the renormalization of the 
coupling. The number of degrees of freedom, Nd. j. in our fits varied between = 344 

and N™* n j = 31 with decreasing f3, from j3 — 10 to ft — 1.015. In all the cases we obtained 
X 2 /Ad.o./.~1-0 (with the only exception being for f3 = 1.015 where x 2 /Nd.o.f. — 2.1). As an 
example we show in Figure || the ^-dependence of the potential Vw(T; R) at T = 1 and 
T = 15 (circles). Crosses show the corresponding values of g^VcouiiT; R) with the fitted 
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Figure 3: The Wilson loop potential Vw(T; R) (circles) and the corresponding perturbative 
potential g|.Vcw(T; R) (crosses) at (3 = 10 on a 16 3 x 32 lattice with and without smearing. 
For T = 15 the potential from smeared and normal Wilson loops would be indistinguishable 
in the figure. 



g\. We would like to emphasize that, even at T = 15, the potential from nonplanar Wilson 
loops cannot be fitted with a continuum-like 1/R ansatz, which can be attempted when only 
planar Wilson loops are available (see e.g. || ITUf). 

At large values of (3 ((3 > 5) our one-parameter fit works well for the whole set of 
data points, i.e. all points with 1 < T < 15 and 1 < R < 8-^/3 (N do j^ = 344). The 
non-Coulomb-like corrections are therefore negligible despite our rather small statistical 
errors. With decreasing (3 these corrections become noticeable at small values of R (recall 
the fast fall-off of V n coui like 1/-R 3 ) and small values of T. After excluding the data points 
corresponding to these small R and T values from the fit, we obtained fits with a high 
confidence level, but fewer degrees of freedom, as described above. Remarkably, the lattice 
Coulomb potential works, at least at large distance, really well even at /3-values very close 
to the phase transition point, (3*, to the confined phase. This observation in fact constitutes 
one of the main results of this paper. 

The extracted values of the renormalized coupling g\ as well as the perturbative values 
9i-ioo P an d gi-ioop are listed in Table [I] and shown in Figure (|. The difference between g\ 
and the perturbative value g\_i oov becomes noticeable at (3^2.0. This is the region where 
one expects the perturbative three-loop contributions to become important. In addition 
nonperturbative effects may start to significantly contribute. 
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In numerical simulations of QCD it has become standard to use a smearing procedure to 
improve the signal-to-noise ratio. In the computation of Wilson loops the space-like link 
elements Uk are replaced be recursively computed "smeared links" 



U 



(n+l) 



xk 



P 



all. 



(n) 



xk 



(n) TT {n) 



ttV) Tj( n )\ _|_ rrWT jjK n ) jjV n ) 
x] x+j,k : x+k,j ' x—j,j u x—j,k x—j+k,j) 



r(»l)t 7-7-H 



r(«) 



j,* = 1,2, 3, (25) 



where P indicates a projection onto the gauge group. The choice of a, the weight of the direct 
link relative to the "staples", and the number of smearing iterations n sme are parameters 
that can be optimized. 

At lowest perturbative order the recursive relation eq. (^3p becomes for Fourier trans- 
formed fields 



.4 



(n+l) 



(P) 



1 



a + 4 



Ar (P) (a + E 2 cos ^) + 4 sin f E (P) 



(26) 



This relation can easily be iterated and inserted in the numerical computation of the 
momentum sums to obtain the lowest order perturbative contribution to the potential from 
the space-space parts of the Wilson loops. 

In Figure ^ we show the effect of smearing on the potential extracted at T = 1 (circles) for 
13 — 10. The parameters of smearing were chosen as a = 4 and n sme = 6. The crosses on the 
top of these circles indicate the values of the potential obtained from the one-parameter fit 
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p 


n 2 

' ' 1 b VJKJIJ 


Qj—Jnnn 


9r 


1.015 


1.22789 


1.33747 


1.892(5) 


1.020 


1.22068 


1.32866 


1.815(5) 


1.030 


1.20652 


1.31138 


1.716(4) 


1.050 


1.17914 


1.27812 


1.585(2) 


1.100 


1.11570 


1.20179 


1.3913(10) 


1.200 


1.00694 


1.07325 


1.1707(7) 


1.500 


0.77778 


0.81173 


0.8379(5) 


2.000 


0.56250 


0.57682 


0.5833(3) 


5.000 


0.21000 


0.21092 


0.21104(4) 


7.000 


0.14796 


0.14829 


0.14834(3) 


10.00 


0.10250 


0.10262 


0.10263(3) 



Table 1: The renormalized coupling g\ extracted from the fit on a 32 x 16 3 lattice as compared 
with 1-loop and 2-loop perturbative values. 

to the perturbative form, as described above. In this case the x 2 ~ va lue is somewhat larger 
(x 2 /Nd.o.f. ~ 2.9). However the extracted value of the renormalized coupling g\ appears 
to be very stable with respect to the smearing procedure, and g\ from the smeared loops 
coincides (within errorbars) with the coupling found without smearing. 

For non-abelian lattice gauge theory it has proved useful, in perturbative computations, 
to use an the 'improved' coupling, such as g 2 E JTJ, [T^| 

g% = -(l-(U p )) , (27) 

where U p denotes the plaquette and c\ = 1/4 is the first coefficient in its perturbative 
expansion. For pure gauge U(l) the perturbative expansion of (U p ) is known to three 
loops00 

(U p ) = l-^g 2 - ^ - 0.0131185/ - 0.00752/ + . . . . (28) 

From this we get a perturbative expansion of g% in terms of the bare g 2 . We then can express 
g 2 in terms of g\ and substitute in the perturbative formula for g\_i OCfp -, eq. (p~6|) , 

9 2 2-ioo P ,E = 9l(l + \gl + 0.0308594y| + ...). (29) 

The plaquette averages, g 2 E and the corresponding prediction for g\_\ OOVtE from eq. ( p9|) is 
listed in Table 0. Comparing with Table |I| we see, that g^-ioopE based on the 'improved' 
coupling g 2 E is indeed a somewhat better prediction for g\ than the purely perturbative 
92-ioo P - This fact can also be seen in Figure |. 

4 We have checked the two-loop coefficient (i.e. ~ g 6 ) which we use below 
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p 


(Up) 


9e 


Qj—lnnn K 


1.015 


0.66147(2) 


1.3541 


1.6800 


1.02 


0.66749(3) 


1.3301 


1.6429 


1.03 


0.67680(2) 


1.2928 


1.5859 


1.05 


0.69104(2) 


1.2358 


1.5002 


1.1 


0.71674(2) 


1.1330 


1.3501 


1.2 


0.751625(10) 


0.9935 


1.1551 


1.5 


0.812599(6) 


0.7496 


0.8363 


2 


0.864810(8) 


0.5408 


0.5835 


5 


0.9486330(24) 


0.20547 


0.21109 


7 


0.9636091(13) 


0.14556 


0.14833 


10 


0.9746754(8) 


0.10130 


0.10262 



Table 2: The average plaquette, (U p ), the extracted g% and the corresponding prediction for 
9l-ioop,E fr o m eq. © 



4 Universal maximal coupling? 

In paper 
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the existence of a universal finite value 47ra c = g c was predicted at the 
deconfinement point, /3*, based on an analogy with the 2d Kosterlitz-Thouless transition. 
The mechanism of the transition in the pure gauge U(l) theory and in the 2d XY model 
was shown to be the termination of the massless phase by a topological disordering. The 
approach to this maximum finite coupling was conjectured as 



gl-c{(3-(3*) 



A 



(30) 



with the estimates g\ = 4vra c = 1.90 ± 0.10 and A = 0.5(1) |3|. 

We have attempted fits of the form eq. ( |3~0|) to our data. A four-parameter fit to the 
5 data points with the smallest (3 gives (3* = 1.007(3), g\ = 2.4(3) and A = 0.27(8) with a 
X 2 = 0.15 for 1 degree of freedom. Including more data points decreases the estimate of (5* 
and increases the estimate of g 2 c . The above estimate for j3* is somewhat smaller then the 
value f3* = 1.011(2) found in |5|. Constraining j3* to 1.011 a three-parameter fit gives now 
g\ = 2.12(2) and A = 0.38(2) with a \ 2 = 1.20 for 2 degrees of freedom, or, using only 4 
data points, g\ = 2.08(4) and A = 0.43(5) with a x 2 — 0.01 for 1 degree of freedom. Varying 
j3* from 1.009 to 1.013 gives variations somewhat larger than the statistical errors. Our final 
estimates are then 



gl 



2.08 ±0.14 



and 



A = 0.43 ±0.10 



(31) 



with the errors dominated by systematic uncertainties. Obviously, more data points in the 
"critical" region would be desirable to make the systematic errors smaller. 

The existence of g\ as maximal coupling depends on the presence of the monopole induced 
phase transition to a confining phase. When the monopoles are suppressed, this phase 
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transition disappears. This leads to the question of whether in that case larger renormalized 
couplings can be achieved than the above g 2 c . This seems indeed the case as found in [17 



in the context of the U(l) Higgs model with suppressed monopoles. We have repeated 
the measurements of the renormalized coupling for pure U(l) with completely suppressed 
monopoles on an 8 3 x 16 lattice. The results are listed in Table || g\ was extracted exactly 
as for the model with standard Wilson action described in the previous section. Since our 
analysis carefully includes finite size (and finite T) effects, a lattice of size 8 3 x 16 is sufficient 
for our purposes here. In Table ^ we list, for comparison, also renormalized couplings for 
the standard model obtained from an 8 3 x 16 lattice. Comparison with Table [l] shows that 
the residual finite size effects are indeed very small, even close to the phase transition. 



ft 


si(NoMo) 


g 2 R (stand) 


0.0001 


5.19(4) 




0.25 


3.97(2) 




0.5 


2.911(11) 




1.0 


1.445(5) 




1.015 




1.86(3) 


1.02 


1.406(3) 


1.79(3) 


1.05 


1.349(3) 


1.579(8) 


1.2 


1.125(2) 


1.172(2) 


1.5 


0.8315(8) 


0.8362(13) 


2.0 


0.5824(5) 


0.5831(2) 



Table 3: The coupling g 2 R from an 8 3 x 16 lattice for the standard U(l) model, and for the 
model with complete monopole suppression. 



From Table ^ we see that at larger ft (ft > 2.0) monopole suppression has no effect, 
as expected, since monopoles are very much suppressed in the standard U(l) model there. 
For smaller ft, due to the effect from monopoles, g 2 R raises faster for the standard U(l) 
model until reaching its maximal value g 2 at the phase transition point. With monopoles 
completely suppressed, g\ raises more slowly. However, the absence of the confinement phase 
transition, allows g\ to keep growing and to become larger than g 2 . Although the model 
with complete suppression of monopole has no phase transition even at ft — 0, we stopped 
our measurements at a small positive value to ensure that the bare coupling remains real. 

The maximal coupling g 2 , however, should be universal for U(l) actions with a monopole 
induced phase transition, for example the Villain model. The renormalized coupling was 
measured for the Villain model in ||, albeit only on a 6 3 x 4 lattice. We simulated an 8 3 x 16 
lattice and extracted g\ the same way as for the Wilson action. Our results are listed in 
Table |. 

For the Villain action, there is no perturbative renormalization of the bare coupling. All 
the renormalization comes from the monopoles. Indeed, up to ft — 1.0 the renormalized and 
bare couplings are almost identical. Only then do the monopoles start to have an appreciable 
influence and g\ becomes larger than the bare coupling and continues to raise until the phase 
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0.645 


2.027(30) 


0.65 


1.915(24) 


0.655 


1.868(12) 


0.66 


1.816(9) 


0.67 


1.732(6) 


0.68 


1.672(5) 


0.7 


1.575(6) 


0.75 


1.402(2) 


0.8 


1.285(2) 


1.0 


1.003(2) 



Table 4: The coupling g R from an 8 3 x 16 lattice for the U(l) model with Villain action. 

transition to the confined phase is reached. We made a fit of our data to the form eq. (|3"0"D. 
The fit works very nicely for < 0.75. Including all 8 data points gives = 1.23 for 4 
degrees of freedom with 

g 2 c = 2.12(4) , 0* = 0.6457(8) and A = 0.43(4) (32) 

in nice agreement with the results for the Wilson action, eq. (|3lD and confirming the 
universality of g c . The result eq. (^2j) is stable when omitting data points from the fits. 

5 Conclusions 

We have made an analytical and numerical study of compact lattice pure gauge U(l) 
theory in the Coulomb phase. The main point of interest was the study of the non-Coulomb 
contributions to the heavy charge potential. For this purpose we calculated perturbatively 
the heavy charge potential Vp(R) defined from the Polyakov loop correlations in the 2- 
loop (i.e., ~ g 6 ) approximation. This calculation shows that the non-Coulomb contribution 
Vp Coul (R) is much smaller than the Coulomb contribution Vp oul (R) , at least at large 
distances. 

The conclusions obtained within perturbation theory were confirmed by numerical 
calculations of the potential Vw{T; R) defined from planar and nonplanar time-like Wilson 
loops. We used as a fit formula for Vw(T;R) the lowest order perturbative expression 
with the only fit parameter being the renormalized coupling constant g 2 R . These fits worked 
well for all distances at large 0, and for sufficiently big T and R even down to very close 
to the phase transition to the strong coupling confined phase. It is worth noting that it is 
impossible to obtain a good \ 2 using for the fit a continuum potential ~ 1/R. The values 
of the renormalized coupling g\ are stable with respect to a smearing procedure. 

We confirmed, with better accuracy, the conjectured existence of a universal maximal 
coupling for U(l) models with a monopole induced confining phase transition. In the model 
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with complete suppression of monopoles, however, even larger renormalized couplings can 
be reached. 

Our main conclusion is that compact pure gauge U(l) theory can serve equally well as a 
non-compact version to describe the physics of free photons in the weak coupling region. The 
only difference is a finite, i?-independent renormalization of the coupling constant g 2 — ► g\ 
in the compact theory. 
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